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The self-consistent field (SCF) approach to the thermodynamics of dense polymer liquids is based on the 
idea that short-range correlations in a polymer liquid are almost independent of how monomers are connected 
into polymers over large scales. Some limits of this idea are explored in the context of a perturbation theory 
for mixtures of two structurally identical polymer species, A and B, in which the AB pair interaction differs 
slightly from the AA and BB interaction, and this difference is controlled by a parameter a. An expansion 
of the free energy to first order in a yields an excess free energy of mixing of the form az(N)<j>A<t>B, in 
both lattice and continuum models, where z(N) is a measure of the number of inter-molecular near neighbors 
of each monomer in a one-component (a — 0) reference liquid with chains of length N. This quantity de- 
creases slightly with increasing N because the self-concentration of monomers from the same chain around 
each monomer is slightly higher for longer chains, creating a deeper inter-molecular correlation hole for longer 
chains. We present a theoretical analysis of the resulting TV-dependence of local liquid structure, and predict 
that z(N) = z°°[l + /3JV -1 ' 2 ], where N is an invariant degree of polymerization, and (3 = is a 

universal coefficient. This and related predictions are confirmed by comparison to simulations of a continuum 
bead-spring model, as well as to previously published lattice Monte Carlo simulations. We propose a way to 
estimate the effective interaction parameter appropriate for comparisons of simulation data to SCF theory and to 
coarse-grained theories of corrections to SCF theory, which is based on an extrapolation of coefficients in this 
perturbation theory to the limit TV — > oo. We show that a renormalized one-loop theory contains a quantitatively 
correct description of the iV-dependence of local liquid structure that we study here. 



I. INTRODUCTION 

In its most general form, 1 the Flory-Huggins theory of 
polymer mixtures postulates that the free energy of mixing 
per monomer A/ in a binary polymer mixture can be ex- 
pressed as a sum of the form 

A/ = fcT^-^ln0 l + A/ int (0,T) , (1) 

i 2 

where Ni is the degree of polymerization of species i, for 
i = A or B, 4>a — 4> ar, d &B = 1 — <t>A are volume frac- 
tions, kT is thermal energy, and A/; n t(</>, T) is an interaction 
free energy per monomer. The original Flory-Huggins lattice 
model 2 was more specific, insofar as it predicted that A/; nt 
should be of the regular solution form A/j nt = C4>a4>b, 
where C is a temperature independent constant. In order 
to capture the variety of behaviors observed in real poly- 
mer mixtures, however, it has long been understoo d 1,3 ' 4 ' 5 that 
A/i n t must be allowed to exhibit an essentially arbitrary de- 
pendence on temperature and composition. 

This generalization of the Flory-Huggins theory is very 
flexible, but not infinitely so: As long as A/j nt (4>, T) is as- 
sumed to be independent of the chain lengths Na and Nb, 
the theory makes falsifiable predictions about how the phase 
diagrams for a homologous set of mixtures vary with Na and 
Nb, which are generally quite accurate. The assumption that 
A/^t (4>,T) is independent of chain length is the essential 



physical content of the theory. 

The intuitive basis of this assumption is the idea that 
A/i n t (4>, T) represents a contribution to the free energy that 
is sensitive only to the local structure of a polymer liquid. It is 
also assumed that this contribution is separable from contri- 
butions that arise from the entropy of mixing in homogeneous 
mixtures and, in a generalization to inhomogeneous liquids, 
from the distortion of polymer chains over large length scales. 
The generalized Flory-Huggins theory described above is the 
homogeneous limit of a correspondingly general form of the 
self-consistent field theory (SCFT) of inhomogeneous poly- 
mer liquids. The key assumption in this SCFT 6 is that the in- 
teraction free energy density at any point in the fluid depends 
only upon the temperature and the average monomer con- 
centration very near that point, independent of chain lengths, 
chain architecture, or compositions at distant points. 

The original Flory-Huggins lattice theory combined an im- 
plicit assumption of locality (i.e., that local fluid structure is 
independent of chain length) with a much cruder random mix- 
ing approximation. The original theory introduced a random 
mixing approximation for the free energy of a lattice model 
in which monomers of types i and j on neighboring lattice 
sites interact with a potential energy Vij . In the absence of 
vacancies, this yields 

A/ irlt = az\ &t t4>A4>B , (2) 
where a = [vab — {vaa + ^sb)/2], and where zi att is the 



lattice coordination number, i.e., the number of lattice sites 
neighboring each site. 

The random mixing approximation is known to substan- 
tially overestimate the actual energy of mixing for lattice 
models. In simulations of a simple cubic lattice diluted with a 
modest density of vacancies, Sariban and Binder 7,8 found that 
the energy of mixing was roughly half that predicted by an 
analogous approximation for a lattice with vacancies. In more 
recent lattice Monte Carlo simulations of diblock copoly- 
mer melts on a diluted fee lattice, Matsen and coworkers*^ 
also considered a lattice mean-field (i.e., random mixing) ap- 
proximation for the order-disorder transition of symmetric 
diblocks, and found that it predicts a transition temperature 
more than twice that observed in their simulations, again in- 
dicating a large overestimation of the energy arising from AB 
pair interactions. 



A. Perturbation Theory 

For simple symmetric models, a much more accurate es- 
timate for Af in the one-phase region can be obtained from 
thermodynamic perturbation theory. We focus here on a class 
of structurally symmetric models of A-B polymer blends, on- 
or off-lattice, in which A and B chains are structurally identi- 
cal, but in which the potential energy for AB pair interactions 
is slightly different than that for AA and BB interactions. 
Many of the coarse-grained models used in previous simula- 
tions of polymer blends fit this description, including some 
used in lattice Monte Carlo (MC) 7 ' 8 '"' 12 '" off-lattice MCiii 5 - 
and molecular dynamics (MD) 16 simulations. 

To construct a perturbation theory, let a be a small parame- 
ter, with units of energy, that is proportional to the magnitude 
of the difference between AB and AA interactions. Here, we 
consider an expansion of the free energy as function of a at 
constant temperature, in a liquid containing structurally iden- 
tical A and B chains of equal length N. The limit a = 
is thus a strongly correlated one-component liquid. At the 
critical composition <j> = 1/2, phase separation occurs when 
a exceeds a critical value a c that decreases as 1/N with in- 
creasing N. An expansion of Af for a < a c to first order 
in a is thus expected to become increasingly accurate with 
increasing N. 

In Sec. |III]of this paper, we show that the first order expan- 
sion for the free energy of mixing is given by a function of 
the form 



A/ 



(3) 



for both lattice and continuum models. In this expansion, the 
coefficient z(N) is an "effective coordination number" whose 
value is sensitive to local correlations in the one-component 
reference state. In a lattice model with nearest-neighbor in- 
teractions, z(N) is found to be equal to the average number 
of inter-molecular nearest neighbors, i.e., the number of lat- 
tice sites neighboring a test monomer that are occupied by 



monomers from a different chain than the one containing the 
test monomer. An analogous definition is given for z(N) in a 
continuum model. 

The main source of error in the random mixing model is the 
neglect of the consequences of the so-called correlation hole. 
The immediate environment of any monomer in a dense poly- 
mer liquid is crowded with other monomers from the same 
chain. In a nearly incompressible liquid, this causes a com- 
pensating depression in the number of neighboring monomers 
from other chains, leading to a correlation hole in the inter- 
molecular distribution function, and thus to a decrease in the 
intermolecular interaction energy. The accuracy of a pertur- 
bation theory for the free energy of a system of long chains 
is a result of the fact that the structure of this correlation hole 
at all a < a c is very similar to that in the a = reference 
liquid, because of the smallness of a c . 

Several authors have previously proposed approximations 
that are either equivalent to first order perturbation theory, 
or very closely related. Miiller and Binder 13 proposed a 
"modified Flory-Huggins" approximation for the free en- 
ergy of mixing Af of a simple lattice model that is com- 
pletely equivalent to the first-order perturbation theory of 
Eq. (0. In discussions of the results of continuum bead- 
spring simulations, both Grest and coworkers 16 and Escobedo 
and de Pablo 15 have proposed approximations for structurally 
symmetric continuum models that are either equivalent^ 6 - or 
nearly equivalent 15 to the perturbation theory discussed here. 
None of these authors explicitly described their proposed ap- 
proximations as a form of perturbation theory. 

Miiller and Binder— and Escobedo and de Pablo 15 have 
both shown that these perturbative approximations provide 
extremely accurate predictions for the critical temperature T c 
observed in simulations of symmetric binary blends. In both 
of these studies, the authors showed that critical point found 
from finite-size scaling analysis of simulation results is very 
well approximated by the expression 



az(N)N 
kT c (N) 



(4) 



that is obtained by using the first-order expansion of Eq. 
(f3]) to predict T c . A somewhat stronger statement can be 
made about the lattice MC data of Binder, Mueller, and 
coworkers 12,13 ' 17 : Their results for the dimensionless parame- 
ter az(N)N/kT c (N) deviate systematically from 2 for finite 
chains, but appear to extrapolate in the limit A*" — > oo to an 
asymptote that is equal to 2 to within statistical uncertainties. 
This convergence was not emphasized by Binder and cowork- 
ers, but appears to be a natural consequence of the identifica- 
tion of modified Flory-Huggins theory as first order perturba- 
tion theory, and the fact that deviations from any low order 
perturbation theory for T c must decrease with decreasing a c , 
and thus with increasing N. 
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B. Chain Length Dependence of Local Correlations 

An important difference between the generalized Flory- 
Huggins theory of Eq. (Q} and perturbation theory of Eq. 
©, which are superficially very similar, is that the coefficient 
z(N) in Eq.® actually does depend slightly upon the chain 
length N. Miiller and Binder 13 obtained values for z(N) in 
bond fluctuation model simulations of one-component melts 
for several values of N, ranging from N =20,. . . ,160. They 
found that z(N) approaches a finite asymptote z°° as N — > 
oo, and that deviations of z(N) from z°° are well approxi- 
mated by a function of the form 

z(N) = z°°[l + {3N- 1 / 2 } . (5) 

where N = Nb 6 /v 2 , b is the statistical segment length, v is 
the volume per monomer, and (3 is a constant that Miiller and 
Binder determined empirically. 

In Sec. |lV]of this paper, we analyze the chain length de- 
pendence of the short distance behavior of the intermolecular 
radial distribution function, and of related quantities such as 
z(N). Our analysis starts from the (verifiable) assumption 
that, in a nearly incompressible liquid, the total radial dis- 
tribution function, including both intra- and inter-molecular 
contributions, changes extremely little with changes in chain 
length N. The intramolecular correlation function does, how- 
ever, change slightly with N, due to changes in the number of 
chemically distant monomers from the same chain that con- 
tribute to the self-concentration around any test monomer. 
This causes a systematic decrease in the depth of the inter- 
molecular correlation function with decreasing N, and thus 
an increase in z(N), as as chemically distant intramolecu- 
lar neighbors are simply replaced by intermolecular neigh- 
bors in the immediate environment of any test monomer. The 
relevant concentration of chemically distant intramolecular 
neighbors of any test monomer is calculated using a simple 
random walk model. 

We predict that, while z°° is a model dependent quantity, 
the coefficient (3 in Eq. <(5j has a universal value 

/g\ 3/2 

/?=(-) =2-64 (6) 

for any structurally symmetric model. 

In Sees. IVland lVTl we test our assumptions by comparing 
the resulting predictions to simulation results. In Sec. |Vj we 
verify the accuracy of Eq. © by comparing it to the lattice 
MC results of Miiller and Binder. In IVI1 we present a more 
extensive comparison to the results of our own off-lattice sim- 
ulations, which allow us to more directly test the assumptions 
underlying our analysis. 

C. Comparing Blend Simulations to SCFT 

Simulations of dense polymer liquids are providing in- 
creasingly precise tests of the assumptions underlying SCF 



theory. Lattice Monte Carlo and continuum simulations of 
simple coarse-grained models have been used to quantify 
slight deviations from the random walk model for polymer 
statistics in melts, and deviations from the RPA description of 
composition fluctuations in both polymer blends 7,12 ' 13 ' 16 and 
block copolymer melt a 10 ' 18 ' 19 ' 20 ' 21 . 

In order to compare either simulation or experimental data 
to SCF predictions for multicomponent systems, however, 
one must somehow choose values for the SCF interaction free 
energy A/j nt ((/), T) and/or the effective interaction parameter 

_ -1 dAfto&T) 
Xe 2kT dej) 2 y ' 

relevant to the RPA analysis of long-wavelength scattering. 1 
[Note that, throughout this paper, / and A/ denote free en- 
ergies per monomer, though we use the same symbols else- 
where to denote free energies per volume. 6 .] Because the 
relationship between A/; nt ((/>, T) and the underlying micro- 
scopic parameters is never known a priori, the temperature 
and (sometimes) composition dependence of A/j nt or \e 
have thus far been determined by fitting RPA predictions to 
the available measurements of composition fluctuations, in 
the analysis of either experiment or simulations. The uncer- 
tainty introduced by this fitting procedure becomes a poten- 
tially serious problem, however, when one's goal is to pre- 
cisely quantify small deviations from RPA predictions, which 
is necessary in order to test theories that predict corrections 
to SCFT, such as the one loop theory. For this purpose, it is 
would be very useful to have independent way of unambigu- 
ously defining and accurately calculating A/i nt for a simula- 
tion model, using the microscopic information that is avail- 
able in a simulation. Here, we propose a way of doing this 
for symmetric models, which is based on an extrapolation of 
perturbation theory to the infinite chain limit. 

SCFT has long been believed to be, in some sense, ex- 
act in the limit of infinitely long chains. The basis for 
this belief is, in part, the predictions of a variety of closely 
related one-loop theories of fluctuation effects. The one- 
loop theor y 6,22 ' 23 ' 24,25,26 is a coarse-grained theory, that, when 
properly interpreted 6,26 , predicts small corrections to the free 
energy of an underlying SCF theory. The relative magnitude 
of the fluctuation correction to the free energy is found 6 ,23 ' 24 ' 26 
to decrease as N^ 1 / 2 with increasing chain length N, imply- 
ing that SCFT becomes exact as N — > oo. This idea makes 
sense, however, only if it is understood that the definition of 
SCF theory used in the construction of the renormalized one- 
loop theory 6,26 is one in which A/j„t is defined by a process 
of extrapolation to the infinite chain limit: That is, it makes 
sense to say that the true free energy / convergences to some 
form of SCF free energy as N — > oo only if the relevant SCF 
is understood to be one whose parameters reflect the liquid 
structure of a hypothetical system of infinitely long chains. 

We propose here that a useful approximation for this 
asymptotic SCF free energy may be obtained for simple sym- 
metric models by considering perturbative expansions of both 
A/int and the true free energy / in powers of a. The idea is to 
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identify the coefficients in the expansion of Afi nt by consid- 
ering the limiting behavior as N — > oo of the corresponding 
coefficients in the expansion of the true free energy for blends 
of finite chains. 

Specifically, we identify the first order contribution to 
A/int as the N — > oo limit of the first order contribution 
to Eq. (0 for Af((j>, N). This yields an SCF energy 

A/ irit ~ az°°<t>A<i>B + 0(o?) (8) 

to first order in a. This first order approximation yields 
the same dependence on composition as the original Flory- 
Huggins theory, but reflects the local correlations present in a 
hypothetical reference system of infinitely long chains. 

In Sec. IVHI we discuss the relationship between the 
O^N^ 1 / 2 ) contributions to z(N) in the first order perturba- 
tion theory for finite N and the corrections to SCF theory 
predicted by a renormalized one-loop theory. To test the com- 
patibility of one-loop predictions with our analysis of the per- 
turbation theory, we consider the predictions of the one-loop 
theory for the value of the derivative d A/ /da in a symmetric 
blend with a = 0. We find that the one-loop prediction for 
this quantity is actually identical to that given by Eqs. (01, ©, 
and Eq. ((5), with the same value for the numerical coefficient 
(3 of the 0(N~ X / 2 ) correction. This implies that the renor- 
malized one-loop theory correctly describes the chain length 
dependence of the correlation hole that is discussed here. 

II. MODELS AND NOTATION 

In what follows, we consider a general class of structurally 
symmetric binary polymer blends, uisng a similar language 
for lattice and continuum models. We consider a system con- 
taining a total of M structurally identical chains, each con- 
taining N monomers, in which 4>aM are of type A and 4>bM 
are of type B. Let a be a small parameter that controls the 
magnitude of the difference between AB and AA interac- 
tions. The state a — is thus a ideal mixture of M physi- 
cally indistinguishable chains, in which a fraction 4>a can be 
chosen to be A chains at random. 



A. Lattice Models 

We consider a class of lattice models in which double occu- 
pancy of lattice sites is forbidden, and in which monomers on 
neighboring sites of types i and j interact with a pair potential 
Vij (a) of the form 

Vij(a) = u + ab i3 . (9) 

Here, u is the interaction between all neighboring monomers, 
a is a small parameter, and fey is a symmetric matrix of 
dimensionless coefficients, with bAB — ^>ba- To main- 
tain the symmetry between the two species, we require that 
bAA — i>BB- The value of the parameter u is relevant if and 



only if the system contains vacancies, because changes in u 
can then effect correlations in the one-component reference 
liquid. 

Binder, Deutsch, and Mulle r 12 i 13 ' 17 have simulated one 
variant of the bond fluctuation model in which (in our nota- 
tion) bAA = and bAB = 1 and another in which bAB = 1/2 
and bAA = — 1/2. 

B. Continuum Models 

We also consider a class of structurally symmetric con- 
tinuum models. Consider a system of M chains of length 
N in a volume V, giving an overall monomer concentration 
c = MN/V or an average monomer volume v = 1/c. The 
total potential energy is the sum of intramolecular bonding 
potentials, which are assumed to be the same for A and B 
chains, plus a sum of nonbonded pair potentials. The pair po- 
tential for monomers of type i and j separated by a distance 
r is assumed to be of the form 

Vij(r) =u{r) + abij(r) , (10) 

with = bBB{r), where a has units of energy. 

In the continum model of Grest et al^, (r) is taken to 
be a purely repulsive shifted Lennard-Jones (LJ) interaction, 
with the same LJ diameter and cutoff distances for all pairs, 
but with an interaction energy that is slightly larger for AB 
pairs than for AA or BB pairs. In this model, &^(r) = 0, 
and u(r) and bAB (r) have the same functional form, differing 
only by a prefactor. A very similar model is used in our own 
continuum simulations, which are presented in Sec. |VI| 

C. Correlation Functions 

To discuss a perturbation theory for continuum models, 
it is useful to introduce some notation for inter- and intra- 
molecular correlations in the the one-component reference 
liquid, with a = 0. Let <?i n ter(r, s, N) denote the intermolec- 
ular radial distribution function (RDF) for a test monomer 
with monomer index s in a reference liquid containing chains 
of length N, defined such that <7i n ter(r, s, N) — > 1 as r — > 
oo. The product c<?i n t or (r, s, N) is thus the probability den- 
sity (probability per volume) of finding any monomer from 
a different chain separated by a vector r from such a test 
monomer. Let w(r,s,N) be an intramolecular correlation 
function in this reference liquid, defined as the probability 
density for finding any other monomer from the same chain 
separated by a vector r from a test monomer with monomer 
index s. Let g to t (r, s, N) be the total RDF for a test monomer 
with a specific monomer index s, defined so that 

cgtot(r, s,N) = u(r, s,N) + cg intcI (r, s,N) , (11) 

and so that gt t(r,s,N) — > 1 as r — > oo. Let 
g intcI (r,N), g tot (r,N), and u)(r,N) denote the averages 
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of ffinter(r, s, N), g to t(r, s, N), and ui(r, s, N), respectively, 
with respect to s. 



III. PERTURBATION THEORY 

We consider a first-order perturbation theory for the free 
energy per monomer 



-kT 
f=^hxZ 



MN 

for either type of model. Here, 



Z = J D[H]e- H ' kT 



(12) 



(13) 



is the partition function, H is the total potential energy and 
J D\R] denotes an integral or (in a lattice model) sum over 
distinguishable microstates. Differentiation of / with respect 
to a yields 



da 



where 



1 /dH\ 
MN \~da~ / 



(14) 



(15) 



To first order in a, the free energy per monomer in a system 
with a composition <f> = 4>a, is thus given by 

kT i 

/(^a)~/ + — ^&ln& + 00(60) (16) 

i 

Here, /o is the free energy per monomer of a corresponding 
one-component reference state, with a — 0, when all of the 
chains are treated as indistinguishable in the calculation of 
entropy. The ideal free energy of mixing term accounts for the 
combinatorial entropy associated with the random labelling 
of chains as A or B. The quantity 9(<fi, a = 0) is evaluated in 
the resulting ideal mixture. 

In the simple case of a lattice model with vaa — u and 
vab = u + a, 9 is simply equal to the total number of AB 
neighbor pairs in the system, divided by the total number 
MN of monomers. 

In a continuum model with bAA = ^>bb = 0, 9 is given by 
the sum of values of bAB (f) for all interacting AB monomer 
pairs, divided by M N. 

The composition dependence of 9 within the ideal mixture 
can be determined by simple combinatorical arguments. To 
show this, we consider lattice and continuum models sepa- 
rately. 



A. Lattice Models 

Let z c {N) be the average number of sites neighboring each 
monomer that are occupied by monomers from a different 



chain, evaluated in the reference state a = (i.e., the av- 
erage number of inter-molecular neighbors per monomer). 
Let w c (N) be the average number of neighboring sites that 
are occupied by monomers from the same chain (the average 
number of intra-molecular neighbors). Let 



y c (N) =z c (N)+w c {N) 



(17) 



be the average total number of occupied neighbors. In the ab- 
sence of vacancies, y c (N) must equal the latice coordination 
number. 

In the reference state, the labelling of chains as type A or 
B is random. The probability that any given pair of inter- 
molecular neighbors will belong to different chains labelled 
with types i and j, respectively, is thus simply <pi<pj. Sim- 
ilarly, the probability that any two intramolecular neighbors 
will belong to a chain of type i is simply <\>i. Thus, for a = 0, 



= lw(N)J24>ihi + lz(N)J2> 



(18) 



The pref actors of 1/2 correct for the double counting of pairs. 
A bit of rearrangement yields 



1 



y(N) + z(N)<l>A<i>i 



(19) 



where we have defined 

y(N) 

z(N) 



Vc{N)b AA 
z c (N){b A B 



bAA) 



(20) 



The free energy per monomer is thus given to first order in a 
by 



/ ~ f + ^ay(N) + Af 



(21) 



where Af is given by Eq. (O, with the definition of z(N) 
discussed above. The resulting first order expansion of Af 
is identical to the modified Flory-Huggins approximation of 
Miiller and Binder 



B. Continuum Models 

In a continuum model, for any value of a, 

= f E/ dr( Cl (r) Cj (0))Mr) , (22) 

ij 

where Ci (r) is the instantaneous concentration of monomers 
of type i and v = V/M N is the volume per monomer. In an 
ideal mixture, with a = 0, random labelling of a fraction cf>A 
of the chains as A and the remainder as B yields 

(ci(r)cy(O)) = S ijC Q(r, N)(t> t +c 2 g intcr (r, N)^, . (23) 



5 



Using expression in Eq. d22l yields an expression for 0(<f>, 0) 
of the same form as that given in Eq. ( fT9] l, and a free energy 
of the form given in Eqs. (fJTJ and (|3j, in which the quanti- 
ties y(N) and z(N) are given for a continuum model by the 
integrals 



y(N) = cj dr g tot (r)b AA (r) 
z(N) = c dr 5 intcr (r)A6(r) 



where 



A6(r) = b AB {r)-b AA {r) 



(24) 



(25) 



In what follows, we will also consider the analogous quanti- 
ties 



y(s,N) = cj dr g tot (r,s,N)b AA (r) 

z(s, N) = c / dr g- mteI (r,s,N)Ab(r) , (26) 



for a test monomer with a known monomer index s, with 1 < 

s < N. 



IV. DEPENDENCE OF LOCAL LIQUID STRUCTURE ON 
CHAIN LENGTH 

In this section, we consider how properties of a one- 
component melt that are sensitive to short range correlations, 
such as z(N), depend upon overall chain length N. Our rea- 
soning applies equally well to lattice and continuum mod- 
els, but we will hereafter adopt a notation appropriate to a 
continuum model. To proceed, we first consider how the in- 
tramolecular correlation function u>(r, s, N) depends upon s 
and N, and then consider how this translates into a corre- 
sponding s- and TV-dependence of Winter (r, s, N). 



A. Intra-molecular distribution 

Each monomer in the one-component reference liq- 
uid is surrounded by a concentration o;(r, s, N) of other 
monomers from the same chain, in addition to a concentra- 
tion cpintoi (r, s, N) of monomers from other chains. Let 



W (r,s,7V) = ^P(r,, 



(27) 



where P(r, s' , s, N) is probability density of finding a spe- 
cific pair of monomers from the same chain, with monomer 
indices s' and s, separated by a vector r. The distribution 
uj(r, s, N) is dominated, for small values of r, by contribu- 
tions from monomers for which \s' — s\ is small. As a result, 
for monomers that are far from either chain end (i.e., far from 



s = 1 and s = N), ui(r, s, N) depends only weakly on chain 
length N and index s. 

In the limit N — > oo, P(r, s, s', N) approaches a function 

P°°(t,\s-s'\)= lim P(r,s,s',N) (28) 

jV— »oo 

that depends only on \s — s'\, and u>(r, s, N) approaches a 
function 

w°°(r)= lim Lj(r,s = N/2,N) (29) 

N— »oo 

that is independent of s, for s not too close to either chain 
end. 

For finite chains, uj{r, s, N) contains a small correction to 
w°°(r), which depends upon both s and N. Let 



Sw(t, s, N) = w°°(r) - u(t, s, N) 



(30) 



where ui°°(r) and u>(r, s, N) are evaluated at the same tem- 
perature and concentration. To estimate 6u> (r, s, N), we as- 
sume that infinite chains and long finite chains have very sim- 
ilar conformational statistics. We thus approximate the differ- 
ence 5ui(r, s, N) for monomers that are not too close to either 
chain end by the contributions to ui°° (r) from monomers with 
s' < and s' > N. In this approximation, 

oo 

5u(r,s,N)~ ^°°(r,As)+ ^°°(r,As) , 

s'=— oo S '=N+1 

(31) 

where As ee |s' — sj. The content of this equation is shown 
schematically in Fig. Q] 

For values of As 3> 1, we may further approximate 
P°°(r, As) by a Gaussian distribution 



Pid(r, As) 



I 2ttAs6 2 



3/2 



exp 



3r 2 



2As6 2 



, (32) 



for a continuous random walk, with As = \s — s' | . The use of 
random walk statistics is justified by the fact that, for s»l 
and | JV — s\ S> 1, |s' — s| ^S> 1 for every term in the rele- 
vant sums. In the same limit, we may also approximate sums 
over s' by integrals to obtain an analytic approximation for 
&Jid(r, s, N). Here and hereafter, a subscript id (for "ideal") 
is used to indicate approximations obtained using this ideal- 
ized continuous random-walk chain model. 

The effective coordination number z(s,N) is sensitive 
only to the distribution of monomers that lie within the range 
of the pair potential from a test monomer. To characterize 
how the self-concentrationo;(r, s, N) within this small region 
depends upon s and N, we consider the s- and iV-dependence 
of the distribution <5wid (r = 0, s, N) for a random walk, eval- 
uated at the the position r = of the test monomer. The 
random-walk model yields a deviation 



Su; id (0,s,N) 



( 3 



V27rfe 2 



3/2 



dA 



^3/2 J a 3 / 2 

N-a 



dA 



3 \ 3/2 2 



2tt 



b 3 



,(33) 
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Further corrections that arise from the breakdown of the ran- 
dom walk model near both chain ends are expected to yield 
subdominant contributions of 0(N^ 1 ). 



B. Inter-molecular distribution 

We now consider how the iV-dependence of the in- 
tramolecular correlation function is reflected in intermolec- 
ular correlations. 



1. A simple argument for z(N) 



FIG. 1: Schematic view of Eq. J3U for the difference Scu(r, s, N) 
between the intramolecular function (r) for an infinite chain and 
the corresponding correlation function u(r,s,N) for monomer s 
on a chain of length N. This difference is attributed to the con- 
tributions to u°°(r) of the monomers s < and s > N that are 
"missing" from the finite chain. The concentration of these missing 
monomers near the test monomer (i.e., near r = 0) can be estimated 
using a random walk model, if s is not too near either chain end. As 
the chain length is decreased from oo, chemically distant monomers 
from the same chain are simply replaced by monomers from other 
chains, while leaving the total RDF g to t (r, s, N) almost unchanged. 



where A = \s — s'\. The random walk approximation clearly 
breaks down for s near the chain ends, as expected on phys- 
ical grounds, since Eq. (|33l predicts a 1/ y/s divergence at 
either chain end. 

The free energy will be shown to depend upon on the aver- 
age value 



Using the above random walk model for Soj (r 
yields 



to id (r = 0,iV) 



1 - 

v \ir 



3/2 



1 



ATI/2 



(34) 

0,s,N) 

(35) 



in which N = Nb 6 /v 2 . The quantity v5oj id (r = 0,N) = 
(6/7r) 3 / 2 iV -1 / 2 is the corresponding volume fraction of the 
"missing" monomers in the vicinity of a randomly chosen test 
monomer. 

Note that the integral with respect to s required to calcu- 
late tt>id(r = 0, N) converges, despite the divergence of Eq. 
(l33l for the integrand at both chain ends. This reflects the 
fact that the average with respect to s is dominated not by the 
contributions of a few monomers near the chain ends, but by 
that of many interior monomers. As a result, our use of a ran- 
dom walk model is sufficiently accurate to correctly calculate 
the prefactor of the dominant 0(N~ 1/2 ) to Suj(r = 0, s, AO- 



Each monomer in a dense melt is surrounded by an envi- 
ronment that contains both monomers from other chains and 
chemically distant monomers from the same chain. Near any 
test monomer, the local concentration of chemically distant 
monomers from the same chain is lower than that found for 
an infinite chain by a concentration 8lu(0, s, N), which is cal- 
culated above. We assume that, in an incompressible liquid, 
this decrease in the self-concentration of a finite chain with 
decreasing N must be almost exactly compensated by a cor- 
responding increase in the concentration of inter-molecular 
near neighbors. We assume furthermore that, in a structurally 
symmetric model with a = 0, the contribution of chemically 
distant monomers on the same chain to the overall radial dis- 
tribution function <?tot(r, s, N) within the small region near 
a test monomer that lies within the range of the potential 
is the same, per monomer, as the contribution of monomers 
from other chains. Only intermolecular near neighbors, how- 
ever, contribute to the quantity z(N). We thus expect a de- 
crease in the local self-concentration by a volume fraction 
vSui(0, s, N) is to increase the effective coordination number 
by a fractional amount vSu>(0, s, N), due to a direct replace- 
ment of chemically distant intramolecular neighbors by in- 
termolecular neighbors. This yields a predicted coordination 
number 



z(N) = z°°[l + vSuj{Q,s,N)} 



or 



z(N) 



7T 



3/2 



NV 2 



(36) 



(37) 



where z°° is the limiting value of z(iV) in the limit N — > oo. 



2. A more detailed view 

The same conclusion may be obtained from a set of more 
explicit assumptions about intra- and inter-molecular correla- 
tion functions in a dense liquid. Let g^ t {r) and 3 i ( J ) ter (r) de- 
note the limits of the gtot and <?i n tcr, respectively, in the limit 
N — > oo. The value of z(N) depends only upon the behavior 
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of the intermolecular distribution Winter (r, s, N) for small r, 
for which |r| is less than the range of the pair potential. We 
make the following assumptions about how the small- |r| be- 
havior of g™ t (r, s, N) and Winter (r, s, N) depends upon chain 
length, in a dense melt in which the polymer coil size y/Nb 
is much larger than the range of the pair potential: 

(1) In an almost incompressible liquid, we assume that the 
total RDF <?tot( r , s, N) is almost independent of both s and 
TV, except for monomers very close to one of the chain ends. 
In this problem, this is what is meant when we say that the 
melt is effectively "incompressible". If true, this implies that 

fftot (r, a ,iV)~ ff ~(r) , 08) 

for all monomers except a few near the chain ends. For this 
to be true for all N, however, the decrease in the intramolec- 
ular self-concentration u>(r, s, N) with decreasing N must be 
exactly compensated by an increase in cgi n t 0I (r, s, N). This 
implies 

C<?intcr 

(r,s,N) 

— c 5intcr 

(r)+6uj(r, S ,N) , (39) 

where 5uj(r, s, N) is defined by Eq. ( f30b . 

Corrections to assumption (l38l can arise, even for val- 
ues of s that are far from either chain end, from contribu- 
tions to g^ t (r, s i N) due to correlations between an interior 
monomer test monomer and an end monomer. The resulting 
end-effect corrections (discussed in more detail in the subsec- 
tion) are of order 1/N, and are found to be negligible com- 
pared to the 1 / %/iV corrections that are obtained from the N 
dependence of 6u>(r, s, N). 

(2) The spatial distribution around a test monomer of 
chemically distant monomers from the same chain is assumed 
to closely mimic the local distribution fftnter of monomers 
from other chains, which we can approximate by the RDF 
Sinter ( r ) f° r a s y ste m of infinite chains. We thus assume that 
the distribution 5u> (r, s, N) of the chemically distant "mis- 
ing" monomers (those that must be removed from an infinite 
chain to form the finite chain of interest) is of the form 

5u,(r,s,N)Kg£ teI (r) (40) 

for large N, values of s far from either chain end, and values 
|r| less than the range of the potential. 

(3) The constant of proportionality in Eq. d40b depends on 
an overall concentration of "missing" monomers over a re- 
gion larger than the range of the potential. We assume that 
an average concentration of missing monomers near a test 
monomer can be obtained by using the random-walk model 
for the return probability 5ui(r = 0,s,N). More precisely, 
we assume that 

8u(t, s, N) ~ 5w id (0, s, N)g£ ter (r) (41) 

under the same conditions that Eq. ( 1401 applies. 

By combining assumptions d39l and d4"Tl ). we find that 

«? intcr (r, a, N) ~ Ctcr(r) [ 1 + vSu id (Q, s, N) ] (42) 



Using this approximation yields 

z(s,N)~z°°[l + v6u id (0,s,N)] , (43) 
where Wid(0, s, N) is given by Eq. d33l , where 

z°° = cJdrg£ tet (r)Ab(v) . (44) 

Averaging with respect to N yields Eq. ([37}. 

C. End Effects 

In addition to the dominant C(7V -1 / 2 ) corrections to z(N) 
predicted above, we also expect to find subdominant 0(1/N) 
corrections to both y(N) and z(N) that arise from true end 
effects, due to the perturbation of the liquid structure near 
chain ends. 

First, consider the contribution of end-effects to the total 
correlation function <?tot(r> s), and to the corresponding inte- 
gral y(s, N). Let <7tot( r ; s, s', N) be a distribution function 
for pairs of monomers with specified monomer indices s and 
s', defined so that (c/N)g to t(r, s, s' , N) is the probability per 
unit volume of finding any monomer with index s' separated 
by r from a test monomer with index s, and so that 

fftot(r,s) = — ^ptot(r, s,s',iV) (45) 

s' 

We assume that the deviation of gtot (r, s, s', TV) from g^ t (r) 
is dominated by pairs of values where either s or s' (but not 
both) are near one of the chains, and the other is somewhere 
in the interior of the chain. Reflecting this assumption, we as- 
sume that g tot (r, s,s', N) can be approximated by a function 
of the form 

5 tot(r, s, s', N) = g™ t {r) + Sg cnd (r, s) + Sg end (r, s') (46) 

where 5g end (r, s) is a deviation that is large only for s near 
1 or N, and vanishes for interior monomers. This func- 
tional form assumes that the correction Sg cnd (r, s) that arises 
from s near one of the chain ends is independent of s' 
when s' is an interior monomer, and similarly for the correc- 
tion 5g ond (r, s'). This approximation captures the dominant 
0(N~ 1 ) corrections to gtot(r), but ignores smaller 0(N~ 2 ) 
corrections arising from contributions in which both s and 
s' are near chain ends. Within this 0(N^ 1 ) approximation, 
evaluating the sum with respect to s' yields 

fftot(r, s) = ffto t (r) + Sg cnd (r, s) + —d end (r) (47) 
where 

4„d(r) = ^<5 5cnd (r,s') . (48) 

s f 

Within the above approximation, we thus predict that, for in- 
terior monomers, <7tot( r j s ) should deviate from g^ t (r) by an 
amount that is proportional to 1 /N but independent of s. 
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This approximation for <7tot(r, s) also implies that we ex- 
pect to find 



6 
N 



(49) 



for interior monomers. The 1/N correction to y(s,N) for 
interior monomers is the result of occasional close contact 
between an interior test monomers and end monomers, and 
is assumed to be independent of the monomer index s of the 
interior test monomer. In addition, we expect to see a much 
larger 0(1) deviation from y°° for the last few monomers at 
either chain end. 

Similar reasoning suggests that the quantity z(s.N) for 
s far from either chain end should also exhibit an 0(1/ N) 
contribution, due to close contacts between the interior test 
monomer and end monomers of other chains, in addition 
to the 0(N^ 1 ^ 2 ) correction described above. The same 
reasoning suggests that this 0(1/ N) correction for interior 
monomers should be independent of s, implying a functional 
form 



z(s, N) ~ z°° [ 1 + vSoj id (0, s,N)] + 



in which z°° and 7 are material parameters. 



1_ 

N 



(50) 



V. COMPARISON TO LATTICE SIMULATIONS 

Fig. |2]shows a comparison of theoretical predictions to the 
lattice Monte Carlo results of Miiller and Binder 13 for z(N) 
for two different variants of the bond fluctuation model. In 
both variants of the model, a monomer is taken to occupy 8 
sites within a cubic lattice, from which other monomers are 
excluded. The volume fraction of occupied sites is 50%. The 
top panel of Fig. [2]shows results for z(N) for chains of length 
N =20, . . . , 160 for a model in which each monomer inter- 
acts with monomers that are located at any of 54 neighbors. 
The bottom panel shows results for a model in which each 
monomer interacts only with monomers at the 6 nearest pos- 
sible positions in a cubic lattice. We have compared both sets 
of data to a prediction 



z(N) 



1+/3N 



-1/2 



i/N 



(51) 



in which (3 is given by Eq. ((6), but in which z°° and 7' are 
treated as adjustable parameters. The j'/N term is included 
to account both for the 1 /N contribution to Eq. ( T50b for in- 
terior monomers, and for contributions to the average over 
s arising from 0(1) deviations from z°° for monomers near 
either chain end. 

In the absence of a prediction for the coefficient 0, Miiller 
and Binder fit each of these data sets to a function z(N) = 
z°°[l + flfitN' 1 / 2 ], while treating both z°° and (3f it as ad- 
justable parameters. This yields best fit parameters (3fu = 
2.846 for z 54 and f3 flt = 3.330 for z 6 slightly higher than 



z"(l + $N~" 2 ) + y/N 



z~(l +|3iV ) 



N 



^ 0.38 




FIG. 2: Lattice MC results of Miiller and Binder 13 for z(N) (sym- 
bols) vs. 7V~ 1//2 for two different variants of the bond fluctuation 
model. The top panel shows results for a model in which each site 
interacts with monomers located on any of 54 nearby sites, while 
the bottom shows results for a model in which each site can interact 
with only 6 neighboring sites. In each panel, the solid line is a best 
fit to Eq. J51l >. using the predicted value of /3 = (6/7r) 3//2 . This fit 
yields z°° = 2.109 and 7' = 0.588 for the model with 54 neighbors 
and z°° = 0.3134 and 7' = 0.310 for the model with 6 neighbors. 
Dashed lines show the estimated asymptote z°° [1 + /37V -1 / 2 ], with 
the same value for in order to show the predicted asymptotic 
slope. 



the predicted value of 2.64. The quality of the fit is approxi- 
mately the same with either functional form. 

Our predictions of a univeral value for the asymptotic slope 
in these plots is consistent with this data. Inclusion of the 
1/N end correction is necessary to adequately fit this data for 
modest values of N, however, particularly for the model with 
very short range interactions. 



VI. CONTINUUM SIMULATIONS 

In order to test the predictions of Sec. |lV]in more detail, 
we have also conducted off-lattice Monte Carlo simulations 
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of a simple bead-spring model. 

Our simulations use a model of flexible polymers with 
a short-range repulsive nonbonded pair potential and a har- 
monic bond potential. The non-bonded pair potential u pa i r (r ) 
in our one-component reference liquid is a purely repulsive 
shifted Lennard- Jones potential, of the form 

Upair(r) = eF(r) 

F(r) = 4[(a/r) 12 - [a/rf] + 1 (52) 

for r less than than a cutoff r c = a2 1 ^ 6 , and v p3l i T (r) = 
F(r) — for r > r c . The bond potential is a harmonic 
spring, 

VbondM = -«(r - I) 2 . (53) 

All of our simulations have used parameter e = kT, I = a, 
and k = 400fcT. 

We have simulated chains of of length TV = 16, 32, 64, 128, 
and 256 at a fixed monomer concentration of c = Q.7a~ 3 . For 
these parameters, we obtain an asymptotic statistical segment 
length b = 1.335, which was obtained by extrapolating the 
slope of a plot of mean squared end-to-end vector versus TV to 
N = oo. All simulations use a cubic LxLxL simulation cell 
with periodic boundary conditions. We have simulated sys- 
tems containing M = 1 176, 588 294, 146, and 146 chains for 
TV =16, 32, 64, 128, and 256, respectively. This correspond 
to L ~ 29.955cr for TV < 128, and L ~ 37.7 for TV = 256. 
The Monte Carlo simulations reported here were carried out 
using a combination of hybrid Monte Carlo / Molecular Dy- 
namics (MC/MD), reptation, and double-rebridging 27 moves. 

Fig. [3]shows the distribution functions g to t (r) and 

Winter (r) 

in the one-component liquid for chains of length TV = 
16, ... , 256. Results for gtot(r) for chains of different length 
are indistinguishable at the scale of the main plot, but the cor- 
relation hole in Winter ( r ) becomes visibly deeper with increas- 
ing TV. The slight dependence of <?tot(r, N) on TV is visible 
in the inset. 



A. Intermolecular Coordination Number 

We consider a perturbation theory for a model of blends 
in which the nonbonded pair potential is of the form used by 
Grest and Lacasse 16 . In this model, the pair potential for i and 
j monomers is taken to be a repulsive LJ potential Vij (r) = 
CijF(r), where F(r) is given by Eq. ( 1521 , in which eaa = 
ebb = e and €ab = e + a, while the same LJ diameter a 
and cutoff distance r c = 2 1 / 6 <r is used for all i and j. For this 
model, 

z(s, N) = cjdr 5 i„tcr(r, s, N)F(r) , (54) 

and z(N) is the corresponding average with respect to s. 

Fig. 2]shows our results for z(N). To compare this data to 
theoretical predictions, we have fit values of z(N) for all but 
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FIG. 3: Intermolecular and total radial distribution functions 
ffinter(r) and gtot(r) for chains of length TV =16,. . . , 256. Results 
for g t ot (r) for chains of different length are indinstinguishable in the 
main plot. Inset: gtot (r) in an expanded scale, in which the slight 
dependence on N is visible. 




FIG. 4: Simulation results for z(N) (symbols) vs. TV -1 / 2 , com- 
pared to the prediction of Eq. ( 15 lb . using the predicted value 
of (3 = (6/V) 3 / 2 . A best fit to the four longest chain lengths 
(TV = 32, . . . , 256), shown by the solid line, yields parameters 
z x = 0.2926 and 7' = 0.0361. The dashed line is the asymp- 
totic line 2°°[1 + /3TV" 1/2 ], with the same values for z°° and (3. 
Error bars are shown, but are very small on this scale. 

the shortest chain (TV = 16) to Eq. ( BTT l. We chose to exclude 
the data for TV = 16 from this fit because the predictions are 
an asymptotic expansion that is expected to be acccurate only 
for sufficiently long chains, and because excluding the short- 
est chains from this data substantially improved the quality of 
the fit. The fit for the 4 longest chains, TV =32,. . . ,256, agrees 
with the data to within our (very small) statistical errors. 

Fig. [5] shows a corresponding comparison of theoretical 
predictions to simulation results for the quantity 

3JV/4 

z"» rf (TV)^__ z(s,N) . (55) 

s=iV/4+l 
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(N) (symbols) vs. N com- 

QTnid 



FIG. 5: Simulation results for z" 

pared to the prediction of Eq. l !56t . using Eq. d57b for /3 m ' a . A best 
fit, shown by the solid line, yields parameters 2°° = 0.2937 and 
7 = 0.1014. The dashed line is the asymptote z°°[l + (3N~ 1/2 ]. 



This is the average of z(s, N) over the middle half of each 
chain. This quantity, unlike the average z(N) over all 
monomers, excludes contributions from monomers very near 
the chain ends. By using Eq. d50l > for z(s,N), and ap- 
proximating the sum over s in Eq. (EBl by an integral over 
N/A < s < 3JV/4, we obtain a predicted TV-dependence 



rnid 



(TV) 



1 _|_ pmidjy-l/2 



1_ 

N 



in which 



f3 mid = (V3-1) (^j 



3/2 



1.932 



(56) 



(57) 



The approximation of a sum over s by an integral gives rise to 
errors of 0(N- 3 / 2 ), which lie beyond the 0(1/N) accuracy 
of Eq. (T50b for the summand z(s, N). Because the sum over 
s that defines z mld (N) only includes interior monomers, for 
which we expect Eq. d50l l for z(s, N) to be valid, we may 
identify the constant 7 in this fit with the constant 7 in Eq. 
(IBTfl l. The prediction fits the data for all TV = 16, ... , 256 to 
within the (small) statistical errors. A fit of the same data to 
z°° [l+fiffiN- 1 / 2 ], in which 0$* is treated as an adjustable 
parameter, yields a slightly worse fit, and a value fi mld = 
2.1233 slightly higher than the predicted value of mld . 

Fig. [6] shows our simulation results for z(s,N) 
for monomers s — 1, . . . , N/ 2 on chains of length 
N =16,. . . ,256. Results for each chain length are compared 
to the predictions of Eq. d50l l. shown by dashed lines. Val- 
ues for the two parameters z°° and 7 have been taken from 
the fit of z mld {N) shown in Fig. [5] Because Eq. ( l33T l for 
<5wid(0, s, N) was derived using continuous Gaussian chain 
model, and yields a value that diverges at both chain ends, we 
have taken 



Suj id (0,s,N) 
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/ 63 
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+ 



(58) 



FIG. 6: Simulation results for z(s, N) (symbols) vs. s, for s 

I N/2 and chains of length iV =16,. . . ,256, compared to the 

prediction of Eqs. J50l l and d58l l. shown by dashed lines. Values for 



the two parameters z° 
shown in Fig. [5] 



and 7 were taken from the fit of z mld (N) 



to apply Eq. (|33T > to a discrete chain with monomer indices 
s = 1, . . . , N, with a = 1/2. Different choices for the value 
of the offset a yield predictions that differ by a correction of 
0(N~ 3 / 2 ). The value a = 1/2 was chosen to minimize the 
discretization error arising from our use of integral to approx- 
imate a sum over s in the derivation of Eq. (1551 . 

Agreement between this data for z(s, N) and Eq. (l50t is 
excellent for all TV, and for all s except values near the chain 
ends, for which all analytic approximations are expected to 
fail. 



B. Total Coordination Number 



We have also considered a quantity 



y(s,N)=c / dr g tot (r,s,N)F{r) 



(59) 



that depends on the total RDF g to t ■ For the blend model con- 
sidered here, for which bAA (r), the quantity y(N) that ap- 
pears in the perturbation theory actually vanishes. We have 
nonetheless considered the quantity defined above as a way 
to test our assumptions about the chain length dependence of 

9intei(r,S,N). 

Fig. [7] shows our results for the average 



2 3N/4 

=N/4+l 



(60) 



of the total coordination number y(s, N) over the middle half 
of each chain, and for y(s, N) itself, respectively. Results 
for y mld (N) have been fit to a predicted form y mld (N) — 
y°° + 8/N, which follows immediately from Eq. (|49l for 
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FIG. 7: Simulation results for y mid (N) (symbols) vs. 1/N. A 
best fit to y mid {N) = y°° + 8/N, shown by the solid line, yields 
parameters y°° = 4.0133 and 8 = 0.4872. 
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VII. SCFT AND ONE-LOOP THEORY 

There are good reasons to believe that a generalized form 
of SCFT becomes increasingly accurate with increasing chain 
length, and is exact in the limit of infinitely long chains. The 
strongest theoretical evidence for this hypothesis comes from 
investigation of corrections to SCFT within the context of the 
renormalized one-loop theory 6 , as discussed below. The hy- 
pothesis that SCFT is exact in the limit of long chains is, of 
course, also consistent with the striking success of the Flory- 
Huggins and (particularly) RPA theories in describing exper- 
imental data from mixtures of long finite polymers. 

The renormalized one-loop theory yields a prediction for 
the free energy per monomer of the form 

f = kTY / ^& + fnA<f>,T)+r , (61) 

i 

in which /* is a one-loop correction to generalized Flory- 
Huggins theory. For symmetric models of the type considered 
here, the predicted correction /* is a function of the form 

kT 

/*(<£, Xc AO = ]y#l73/*(6 XeN) (62) 

where /* is a dimensionless function of <fi and XcN. Here, 
Xc is an effective interaction parameter that is related to 
A/i n t(0, T) by Eq. @. The SCF interaction free energy 
/int, which is required as an input to the theory, can have an 
arbitrary composition dependence. To compare one-loop pre- 
dictions to our perturbation theory of symmetric mixtures, we 
must also allow /; nt to depend upon the parameter a of the 
underlying microscopic model. 



FIG. 8: Simulation results for y(s,N) (symbols) vs. s for 
s = 1, . . . , N/2 and N =16,. . . ,256, compared to the prediction 
y(s, N) — y°° + 8/N (dashed lines). Values for the parameters y°° 
and 8 were taken from the fit shown in Fig. [8] 



y(s, N) for interior monomers. Note that the fractional de- 
viations of y mld {N) from y°° are much smaller than those 
found for z mld (N): For the shortest chains, with N = 16, 
y mld (N) deviates from y°° by about 1%, whereas z mld (N) 
deviates from z°° by roughly 30%. 

In Fig. [8] we compare data for y(s,N) for all s and 
N to the functional form predicted in Eq. ( |49l l for interior 
monomers. Here, we have used the values of y°° and 5 ob- 
tained from the fit shown in Fig. Q This data clearly con- 
firms that the small deviations y(s, N) from for interior 
monomers are independent of s and proportional to 1/N, as 
expected if the dominant corrections to y°° arise from oc- 
casional contact of interior monomers with end monomers. 
The inset to Fig. [8]shows the much larger 0(1) deviation of 
y(s, N) from y°° for the last bead at each chain end. 



A. SCFT 

Because the one-loop contribution /* decreases as TV -1 / 2 
with increasing chain length, the SCF contribution f mt may 
be identified by considering the limit N — > oo of the true free 
energy /. As discussed in the introduction, an expansion of 
/int to first order in a may thus be obtained by simply taking 
the limit N — > oo of the corresponding first order expansion 
of /, which is given in Eq.(0l. This yields 



/int =fo+ \ a V°° + azC 



(63) 



An expression for the corresponding contribution Af int to 
the free energy of mixing is given in Eq. (0. This result 
yields an expansion of effective SCF interaction parameter 
Xc, to first order in a, as 



Xc 



kT 



(64) 



The fact that this expansion for Xc is independent of com- 
position is a special feature of the expansion of this class of 



12 



symmetric models to first order in a: We do not expect it to 
survive any generalization to structurally asymmetric models 
or to higher order in a. Nothing rigorous can be said about 
the temperature dependence, even in this first order expan- 
sion, because the expansion with respect to a has been carried 
out at constant T, and z°° thus has an unknown dependence 
on T. The coefficient z°° would be independent of T only in 
an athermal reference system, such as a lattice model with no 
vacancies or a model of tangent hard spheres. 



B. One-Loop Contribution 

Both the one-loop theory for / and the simple perturbative 
expansion of / predict corrections to the SCF free energy, as 
defined above, that are of order A~ 1//2 . The one-loop theory 
is simply not a perturbation theory, since it predicts a cor- 
rection f *(</), XeN) that is a nonlinear function of XeN, and 
that exhibits singular behavior near the spinodal. We can test 
whether one-loop predictions are consistent with our analysis 
of a more microscopic perturbation theory, however, by con- 
sidering the predictions of the one-loop theory for the deriva- 
tive 8 = df I da at a = 0, and comparing expressions for the 
C(iV _1 / 2 ) contribution to this coefficient. 

The one-loop expression for 9 = df /da at a = can be 
expressed as a sum 



Of _ dfint | df* 
da da da 



(65) 



The SCF contribution is simply 



dfi, 



da 



a=0 



~y°°+z°°<j> A <pB 



(66) 



and 9* = df* /da. The one-loop correction is of the form 

df* _ d Xc (<P,a) df*(4,, X e,N) 



da 



da 



dXc 



(67) 



It is straightforward to show, by using the functional form 
given in Eq. d62| i, that this correction is proportional to 

#-1/2. 

In the accompanying paper, we explicitly calculate the re- 
quired derivative of the one-loop correction /*, and find that 



dr_ 

da 



vSu id (r = 0,N)</)A<t>i 



(68) 



a=0 



Combining this with Eqs. ( 1651167b and Eq. d64l ) for Xc yields 
1 



dj_ 

da 



-y°° + z™<j>A<t>B [1 + v8Q id {r = 0, N)] 



This is identical to the expression obtained in Sec. [IV] 

We thus conclude that the one-loop theory implicitly con- 
tains a correct description of the TV-dependence of the cor- 
relation hole. This is enough to guarantee that the one- 
loop theory will yield a very accurate description of correc- 
tions to SCFT in weakly non-ideal symmetric blends, with 
X C N <C 1. In an accompanying paper, we confirm that this is 
true by comparing one-loop predictions to simulation results 
for composition fluctuations in such blends. There, we also 
examine the accuracy of the one loop theory for larger values 
of a, up to the critical value, for which perturbation theory 
becomes inadequate. 



VIII. CONCLUSIONS 



A simple physical picture has been given for how intra- and 
inter-molecular correlation functions vary with chain length 
in a polymer melt. A theory based on this picture is in 
excellent agreement with computer simulation results. The 
structure of a one-component melt is related by perturbation 
theory to the free energy of mixing in corresponding struc- 
turally symmetric blends. The ©(A -1 / 2 ) contribution to the 
depth of the intramolecular correlation hole in the melt of fi- 
nite chains leads to a slightly higher free energy of mixing 
in mixtures of shorter chains. This simply reflects the fact 
that monomers on shorter chains are less strongly screened 
from contact with other chains. Perturbation theory may be 
used to estimate the SCF interaction free energy appropriate 
for comparison of SCF theory to simulations, by identifying 
SCF theory with the A — > oo limit of the perturbation the- 
ory. If this prescription is used identify SCF parameters, the 
predictions of the one-loop theory for corrections to SCF the- 
ory is found to be consistent with the perturbation theory pre- 
sented here, insofar as both theories give identical results for 
a 0(N^ 1 ' 2 ) correction to the apparent interaction parameter 
in weakly non-ideal symmetric mixtures, with %A -C 1. 



(69) 
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